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Abstract. We study the influence of the symmetry en- 
ergy of the equation of state on the thermal evolution of 
protoneutron stars and the properties of supernova neutri- 
nos by the numerical simulations after the protoneutron 
stars are formed. As for the equation of state (EOS) of 
nuclear matter, we take two EOS's with different symme- 
try energies obtained by the relativistic mean field theory. 
We find the symmetry energy plays the essential role on 
the evolution of lepton profiles and the neutrino fluxes. 

Key words: stars: neutron - supernovae: general - dense 
matter - equation of state 



1. Introduction 

The recent advance in radioactive nuclear beam experi- 
ments provides us with novel information on unstable nu- 
clei far away from the stability line (Tanihata et al. 1985). 
The symmetry energy becomes more essential to under- 
stand the nuclear structure as we go further away from 
stability. The symmetry energy is also important to pro- 
vide the equation of state (EOS) for neutron stars and 
supernovae. The properties of dense matter under neu- 
tron rich environment determine the structure and chem- 
ical composition of stars and may change their evolution 
drastically. Therefore, it is very interesting to study the 
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influence of the symmetry energy on astrophysical prob- 
lems. 

The effects of the symmetry energy on the evolution of 
neutron stars and supernovae have been studied by several 
authors. The possibility of the rapid cooling of neutron 
stars due to the direct URCA process has been discussed 
in the case of the large proton fraction of neutron star 
matter, which is sensitive to the density dependence of 
the symmetry energy (Boguta 1981; Lattimer et al. 1991). 
Bruenn (1989) studied the effects on supernova explosions 
systematically by doing numerical simulations of gravita- 
tional core collapse with the parameterized EOS. Swesty 
et al. (1994) also studied the role in the prompt phase of 
the supernova explosion. However, there has been no sys- 
tematic study on the effect of the symmetry energy on the 
birth of neutron stars and supernova neutrinos as far as 
we know. Moreover, in most of previous studies, the pa- 
rameterized formula have been used to provide the EOS of 
dense matter beyond the normal nuclear matter density. 

Recently, there has been a great progress in the 
study of nuclei and dense matter within the relativistic 
many body framework (Serot & Walecka 1986). It was 
demonstrated that the relativistic Brueckner Hartree Fock 
(RBHF) theory is capable of reproducing the saturation 
property of nuclear matter starting from the nucleon- 
nucleon interaction determined by the scattering exper- 
iments (Brockmann & Machleidt 1990). The relativistic 
mean field (RMF) theory has been shown to be very suc- 
cessful as an effective theory to describe the ground state 
properties of nuclei in the wide mass range of the periodic 
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tabic (Gambhir ct al. 1990) and has been applied to the 
EOS for neutron stars (Serot & Walecka 1986). It is amaz- 
ing that the RMF theory describes also the properties of 
unstable nuclei away from stability extremely well (Hirata 
et al. 1991; Sugahara & Toki 1994). Having the framework 
constrained by unstable nuclei, Sumiyoshi and Toki ap- 
plied the same RMF theory to provide the data table of 
the EOS for neutron stars and supernovae (Sumiyoshi & 
Toki 1994; Sumiyoshi et al. 1995), which enables us to do 
numerical simulations of the thermal evolution of neutron 
stars and supernovae quantitatively while taking care of 
the experimental data of unstable nuclei. The properties 
of unstable nuclei have been shown to be very sensitive to 
the symmetry energy in the RMF theory (Sumiyoshi et al. 
1993a). Hence, it would be nice if we could see how sen- 
sitive are the properties of neutron stars to the symmetry 
energy and what is the influence of the symmetry energy 
on the birth of neutron stars and supernova neutrinos. 

The numerical simulations of the cooling of protoneu- 
tron stars and neutrino burst have been done by several 
groups (Burrows & Lattimer 1986; Burrows 1988; Suzuki 
1993). The influence on the numerical simulation of the 
birth of neutron stars due to the difference of EOS was 
studied by adopting the EOS tables in the two different 
many body frameworks (Sumiyoshi et al. 1993b). The ef- 
fects of various thermodynamical properties were pointed 
out there besides the stiffness of EOS, which has been 
mainly studied so far (Burrows 1988). Here, we focus on 
the influence of the symmetry energy by changing solely 
the strength of the isovector interaction in the RMF the- 
ory in order to see its effect clearly. 

The purpose of this paper is to explore the influence 
of the symmetry energy on the cooling of protoneutron 
stars just born in supernova explosions. We construct the 
tables of EOS for supernova simulations in the RMF the- 
ory and make comparisons of EOS's when we change the 
symmetry energy drastically. Then, we perform numerical 
simulations of the birth of neutron stars and supernova 
neutrinos adopting the two EOS's with different symme- 
try energy. We investigate the influence of the symmetry 
energy on the thermal evolution of protoneutron stars and 
the properties of neutrino burst emitted during the cooling 
stage. 

This paper is arranged as follows. In section ^, we 
briefly describe the relativistic EOS for neutron stars and 
supernovae. In section |[ after a short introduction on the 
birth of neutron stars, we present the results of our numer- 
ical simulations. We discuss the effects of the symmetry 
energy on the deleptonization, the thermal evolution and 
the supernova neutrinos in the subsections. We summarize 
this paper in section ||. 

2. Relativistic equation of state 

We start with a brief explanation of the tables of EOS for 
numerical simulations. We calculate all the physical quan- 



tities of dense matter within the relativistic mean field 
theory. We refer the review article by Serot & Walecka 
(1986) as for the relativistic many body framework for 
nuclei and dense matter. All the details on the relativistic 
EOS for neutron stars and supernovae in the RMF theory 
has been reported in the recent papers. (Sumiyoshi & Toki 
1994; Sumiyoshi et al. 1995). 

We adopt the phenomenological lagrangian with the 
non-linear a and lo terms, which is motivated by the re- 
cent success of the RBHF theory (Brockmann & Mach- 
leidt 1990) and has been shown to be very successful both 
for nuclear properties and dense matter (Sugahara & Toki 

1994) . The best parameter set for the lagrangian, named 
TM1, was determined by the least square fitting to a set 
of nuclei including unstable ones, which is important to 
constrain the isovector interaction in the theory. It is re- 
markable that the RMF theory with the parameter set 
TM1 has been demonstrated to reproduce successfully the 
properties of unstable nuclei other than the ones used in 
the fitting. The properties of nuclear matter in the RMF 
theory with TM1 thus constrained has been shown to be 
quite similar to the properties of nuclear matter derived in 
the RBHF theory (Sugahara & Toki 1994). Extending the 
RMF theory to the case at finite temperature, the table 
of the numerical data of physical quantities under various 
conditions of chemical composition, temperature and den- 
sity, which are required for the numerical simulations, was 
constructed for the parameter set TM1 (Sumiyoshi et al. 

1995) . We use this table of EOS as a standard one. 

In order to explore the influence of the symmetry en- 
ergy, we newly construct the table of EOS with a re- 
duced value of the symmetry energy in the RMF theory. 
We reduce the value of the coupling constant g p between 
isovector-vector p meson and nucleon in the lagrangian, 
which is essential to determine the symmetry energy, while 
keeping other parameters of TM1 unchanged. Hereafter, 
we call the modified parameter set as TMS. The symmetry 
energy is a sym = 36.9 MeV for TM1 and a sym = 28.2 MeV 
for TMS: the corresponding coupling constant is g p — 4.63 
for TM1 and g p = 3.50 for TMS. We note that the sym- 
metry energy in TM1 has been checked by unstable nuclei 
and this modification in TMS is to explore the effect of the 
symmetry energy on astrophysical applications. We dis- 
play in Fig. [l] the energy per baryon of symmetric nuclear 
matter and pure neutron matter for the cases of TM1 and 
TMS. The incompressibility K at the normal nuclear mat- 
ter density is 281 MeV for both cases, since the isovector 
meson does not contribute to symmetric nuclear matter. 

We comment here on the density dependence of the 
symmetry energy. The symmetry energy for both cases 
of TM1 and TMS has a monotonically increasing feature 
with density, which is common for all the relativistic many 
body calculations due to the contribution of the isovector 
meson. This feature is demonstrated by the microscopic 
calculation with the use of the RBHF theory (Li et al. 
1992). The reduction of the symmetry energy at the sat- 
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uration density by changing the strength of the isovec- 
tor interaction clearly corresponds to the reduction of the 
symmetry energy at high density. In contrast, the symme- 
try energy in the non-relativistic many body calculations 
such as the one by Wiringa et al. (1988) has generally a 
weak density dependence and has sometimes a decreasing 
feature at high density. A variety of the density depen- 
dence depends on the choice of the density dependent po- 
tential, which is introduced to reproduce the saturation of 
nuclear matter. Therefore, the relation between the sym- 
metry energy at the saturation density and its behavior at 
high density is ambiguous in non-relativistic many body 
calculations. 

When we apply the EOS's in the two cases to neutron 
stars, we found that the chemical compositions are very 
different while the hydrostatic structures are quite simi- 
lar. In Fig. |^-a, we display the proton fraction of neutron 
star matter, which is the ratio between the proton density 
and the baryon density, as a function of the baryon mass 
density. The proton fraction in the case of TMS is smaller 
than the case of TM1 because of the reduced symmetry 
energy. We show in Fig. ^|-b the neutron star mass as a 
function of the central baryon mass density. In contrast 
to the difference in the chemical composition, the neutron 
star masses in the two cases are found very similar. The 
maximum mass turns out to be almost the same value of 
2.2 Mq. The difference in the central baryon mass density 
is only about 1 % for the case of neutron stars with the 
gravitational mass of 1.4 Mq. 

3. Evolution of protoneutron stars 

First of all, we describe briefly the supernova explosions 
and the birth of neutron stars (Suzuki 1994). Gravitational 
collapse of the core of a massive star leads to explosion of 
the envelope (supernova explosion) and formation of the 
neutron star. We shall divide the series of stages starting 
from the onset of the core collapse to the birth of neutron 
star into two phases. 

The first phase is the dynamical phase. The core of a 
massive star becomes unstable when it grows to the point 
near the Chandrasekhar mass (~ 1.4 Mq), and it begins 
to collapse. The collapse never ceases until the central den- 
sity exceeds the nuclear density. Sudden stiffening of the 
EOS above the nuclear density stops the inner core (0.5 - 
0.8 Mq) and the bounce of the inner core launches a shock 
wave into the falling outer core. The matter of the falling 
outer core is swept and decelerated by the shock wave 
and, then, accretes onto the unshocked inner core which 
has been in hydrostatic equilibrium in its dynamical time 
scale of milliseconds. The shock wave expels the envelope 
and we identify it as supernova explosion. The remnant at 
the center which consists of the unshocked inner core and 
the shocked outer core is called as a protoneutron star. It 
contains so many leptons and protons (~ 30 %) that we 



cannot call it a neutron star at this stage. It takes only 
about 1 second for this dynamical process to take place. 

The second phase, which is the quasistatic phase of 
protoneutron star cooling, follows this dynamical phase. 
After the shock wave breaks out of the core surface and the 
accretion onto the inner core ceases, the protoneutron star 
evolves quasistatically keeping hydrostatic configuration. 
Neutrinos which are trapped in the core diffuse out of the 
protoneutron star in a time scale of neutrino diffusion, 
which is of order of 10 seconds. They drive the evolution 
of the protoneutron star into the cold neutron star by 
carrying out thermal energy and lepton number from the 
protoneutron star. Concerning supernova neutrinos, about 
a half of the total energy is emitted during the dynamical 
phase and the rest is emitted during the quasistatic phase. 

In this paper, we study the influence of the nuclear 
symmetry energy on the second stage: the evolution of the 
hot protoneutron star into the normal neutron star. The 
quasistatic cooling of the protoneutron star is simulated 
numerically by solving the general relativistic equations 
for hydrostatic structure (Oppenheimer-Volkoff equation) 
and for neutrino transport with the deleptonization and 
the entropy change of the matter simultaneously using 
the Henyey-type method. As for the neutrino transfer, 
we adopt the multigroup flux limited diffusion scheme 
(Bruenn 1985). We take into account the energy depen- 
dence of neutrino transport coefficients in the multigroup 
scheme. The flux limiter should be introduced in order 
to express the neutrino flux in the transparent regime in 
terms of the diffusion flux. We adopt Mayle and Wilson's 
flux limiter (Mayle et al. 1987) in this work. We include 
the general relativistic effects such as the time dilation and 
the red shift of the neutrino energy. We treat explicitly, v c , 
u e and v^/t, where v^/ T represents the average of v^, P M , 
v T and v r . This is a good approximation in the case where 
we can neglect the existence of \i and r leptons because of 
the low temperature (^ 100 MeV). 

The following neutrino interactions are included as 
opacity sources or collision terms in the neutrino trans- 
fer equations. 



p e < — > n f G , n c + < — > p v c 

A e~ < — > A' v c 

e _ e + * — > v P, plasmon < — ► v D 

N N' < — > N N' v v 

N v < — ► N v, A v < — ► A v 

v < — ► v 



where v represents all species of neutrinos, A is a represen- 
tative heavy nucleus, and N is either a proton or a neutron. 
Most of the interaction rates are taken from Bruenn (1985) 
with some modifications. At present, many body effects on 
neutrino opacity are not included except for the multiple 
scattering suppression effects on nucleon bremsstrahhmg 
process (Raffelt & Seckel 1991). The evolution of the pro- 
toneutron star is driven by the exchange of energy and 
lepton number between the matter and neutrinos due to 
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the above interactions and the neutrino transport. De- 
scriptions of our numerical code are also given in Suzuki 
(1993). 

We perform numerical simulations of the protoneutron 
star evolution using the tables of EOS in the RMF the- 
ory with TM1 and TMS. As for the EOS of the low den- 
sity matter (p < 10 14 g/cm 3 ), Wolff's EOS (Hillebrandt 
& Wolff 1985) is used for both cases. We construct the 
initial models for our numerical simulation by referring to 
the numerical results of Mayle and Wilson at 0.4 sec after 
the core bounce in the hydrodynamical simulations of the 
supernova explosion (Mayle & Wilson 1989; Wilson 1990). 
Initial models for the simulations using the two EOS's are 
constructed using the same entropy profile and the same 
electron fraction profile in order to study the effects of the 
difference in the symmetry energy. Of course, in princi- 
ple, these profiles should differ from each other because 
the difference of EOS should also affect on the dynami- 
cal phase. To extract the direct influence of the nuclear 
symmetry energy on the protoneutron star evolution, we 
neglect its influence on the dynamical phase. 

The calculated initial hydrostatic structures of the 
two protoneutron stars are found similar; the differences 
in the initial densities and the temperature profiles are 
smaller than 3 %. On the contrary, there is a large dif- 
ference in p eq = p p + p c — p n which is the chemical 
potential of v c in the /3-equilibrium with the matter be- 
cause p p and p n are directly affected by the symmetry 
energy. In Figs. we present the profiles of the two 
initial models. p eq at the center of the initial models are 
142 MeV (a sym = 36.9 MeV) for TM1 and 170 MeV for 



TMS fa. 



28.2MeV), respectively. The resultant lep- 



ton fractions, Yl = (n c - + n„ c — n c + — n Po )/nb a 
0.330 for TMS and 0.315 for TM1. 



3.1. Deleptonization 

Starting from these initial models we simulate the evolu- 
tion of the protoneutron stars for the following 15 seconds 
with the same numerical code. Figs. ||-|| also show the 
profiles of the protoneutron stars at the end of the calcu- 
lation (t = 15 sec) for the two models. While the density 
profiles arc still nearly identical, the distributions of the 
temperature and the electron fraction differ largely from 
each other. Especially, the lepton fraction at the center for 
TM1 (0.179) is larger than that for TMS (0.155), while the 
former is smaller than the latter at the initial stage. This 
means that the deleptonization proceeds faster for TMS 
than for TM1. We find that this is caused by the differ- 
ence of p eq . In the central region of protoneutron stars, 
electron type neutrinos are degenerate and their diffusion 
fluxes are roughly proportional to —\ v dn Ve /dr, where A„ 
is the mean free path of v c and n Ue is the number density 
of v e . A„ is roughly inversely proportional to the neutrino 
energy squared, A^ oc (pf, q )~ 2 , and n Ve oc (pf, q ) 3 in the 
degenerate limit. Consequently, the diffusion flux in the 



degenerate limit is proportional to —dp eq jdr; the smaller 
symmetry energy results in the larger p eq , the larger v c 
flux, and therefore the faster deleptonization. 

3.2. Thermal evolution 

The influence of the symmetry energy on the evolution of 
the temperature is complicated in our models. The cen- 
tral temperature at t = 15 sec are 23.2 MeV for TM1 and 
26.3 MeV for TMS. Higher temperature for TMS is due to 
the higher total (matter + neutrinos) entropy per baryon 
(1.29 for TMS and 1.20 for TM1) and the lower lepton 
fraction (0.155 for TMS and 0.179 for TM1). The lower 
lepton fraction which is the consequence of fast delep- 
tonization means the smaller lepton number density and 
the temperature corresponding to a given total entropy 
becomes higher. 

As for the difference in the central entropy, we ana- 
lyze the results of numerical simulations in detail and find 
the following reason for the present case. In the first 10 
seconds of the protoneutron star cooling, the matter en- 
tropy of the central region increases. This entropy change 
has three main origins; inward flux of v^/ T , the down- 
scattering of v c and the emission of v e due to electron 
capture. Since, during the first 10 seconds, the tempera- 
ture profile has its peak in the middle region of the pro- 
toneutron star, not at the center, there is the negative 
gradient in the number density of v^/ T at the central re- 
gion. The flux of v^/t is proportional to the gradient. 
Therefore, at the early stage, v^/t flow inwards in the 
central region and transport the heat into the central re- 
gion from the middle region. The pair annihilations of v^/ T 
(y^v^j 'v T D T — > e _ e + ) increase the matter entropy. In addi- 
tion, in the central region where the electrons are strongly 
degenerate, z/ c -electron scattering (v c e~ — > v c e~) leads to 
the increase of the matter entropy. Neutrinos lose their en- 
ergy at the scattering (downscattering) because the scat- 
tered electrons which were within the Fermi sea should 
have energy greater than the Fermi energy. On the other 
hand, emission of v c (pc~ — > i/ e n) decreases the matter 
entropy. The two heating processes and the one cooling 
process result in the net heating of the central region. 

Detailed analysis of the numerical simulations reveals 
that, among the above three processes which alter the cen- 
tral entropy, the largest difference of heating/cooling rate 
(dS ma t/dt(r — 0)) due to each process between the two 
models is the difference of the cooling rate due to the v e 
emission. Furthermore it is found that the difference in 
dS ma t/dt(r — 0) due to the v c emission is caused mainly 
by the difference of p eq as the case of the deleptonization 
rate. dS ma t/dt due to pe - «-> z/ c n can be expressed as 

dS ma t 



dt 



pe — <->-i/ e n 
1 / du y 



dn v 



pe *-^ e n 



dt 



(i) 



pe <-*i/ e n j 
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where u Ve is the energy density of v c . Since, in the cen- 
tral region at the early stage where the electron capture 
proceeds, dn Vt ,/dt\ pe -<_+ Ven is positive, it can be seen that 
the larger \i e ^ for the smaller symmetry energy results in 
the larger dS ma t/dt due to u c emission, that is, the less 
entropy loss of the matter because dS ma t/ 'dt\ pc - is 
negative there. This small cooling rate due to the v c emis- 
sion leads to the higher entropy and higher temperature 
for the case of smaller symmetry energy. We note, however, 
that the temperature profile is affected by many factors. 

3.3. Supernova neutrinos 

The difference in the evolutions of the protoneutron star 
due to the symmetry energy is reflected with the prop- 
erties of supernova neutrinos emitted during the cooling 
stage of the star. We show in Fig. ^ the net flux of the 
electron type lepton number from the protoneutron star; 
the flux of v e minus u e . The net flux is larger for TMS 
than for TM1 during 15 seconds because of the difference 



in the diffusion fluxes as we have discussed in section 3.1 



This fact corresponds to the higher deleptonization rate 
in TMS having the smaller symmetry energy. Since the fi- 
nal electron fraction of the neutron star at the end of the 
deleptonization is smaller for TMS, the total net neutrino 
number emitted during the cooling is larger for TMS than 
for TMI. It may be possible to extract the information 
of the final electron fraction of the neutron star from this 
quantity. 

We show in Fig. [l(] the calculated time profile of the 
mean energy of v c . The mean energy for TMS becomes 
higher than that for TMI at a later stage. This is due 
to the difference of the temperature profile of the pro- 
toneutron stars, as shown in Fig. ^|. The mean energy of 
neutrinos depends mainly on the temperature at the neu- 
trinosphere. The difference of the mean energies is larger 
for zv/r an d smaller for v c depending on the position of its 
neutrinosphere. The neutrinosphere for u c which interact 
most strongly with the matter locates in the outermost 
region, where both the temperature and the density are 
low and the difference between the temperature in the two 
cases is small. 

We show in Fig. [ll] the luminosity of V Cl which is higher 
for TMS than for TMI. The general feature is quite similar 
for the other types of neutrinos. These results correspond 
to the larger flux and the higher mean energy of neutrinos. 
The total energy carried out by neutrinos is larger for TMS 
than for TMI, since the gravitational mass of the neutron 
star at zero temperature is smaller for TMS. We note that 
the total baryon mass of the protoneutron stars is fixed 
to be 1.62 Mq (the gravitational mass of the initial pro- 
toneutron star is 1.5479 M for TMS and 1.5482 Mq for 
TMI) in the present study and the gravitational mass of 
the cold neutron star turns out to be 1.475 Mq for TMS 
and 1.485 M for TMI. Therefore, the total energy re- 
leased from the protoneutron star amounts to 0.073 Mq 



(1.31 • 10 53 erg) for TMS and 0.063M Q (1.13 • 10 53 erg) 
for TMI, which means the larger luminosity for TMS, 
provided that the time duration of neutrino emission is 
similar. 

4. Summary and Discussions 

The main influence of the symmetry energy of the EOS 
on the birth of neutron stars comes from the change of the 
chemical composition rather than the stiffness. A change 
in the chemical potential for neutrinos affects the rates of 
the interaction with the matter even if the matter density 
is unchanged. It affects the deleptonization rate through 
the change of the diffusion flux, and the temperature in- 
side stars through the change of the heating rate. Those 
changes influence the properties of the neutrino burst such 
as the time profile of the mean energy while the density 
profile remains unchanged. The difference in the proper- 
ties of the neutrino burst is emphasized moreover by the 
difference in the proton fraction and the binding energy 
of the cold neutron stars at the end of the birth stage. 

The above results come from our study on the evolu- 
tion of protoneutron stars by assuming the same initial 
conditions following the result of Mayle & Wilson. Here 
we comment on expected difference in the initial protoneu- 
tron star configuration when we simulate also the dynam- 
ical phase, the birth of the protoneutron stars, using the 
two EOS's. Since TMS has a smaller coupling constant 
<? p , it is expected that the symmetry energy and hence 
the difference between the chemical potential of neutrons 
and that of protons are smaller for TMS than for TMI 
even at the density of the collapsing core. This means 
a larger fraction of free protons for TMS, which will re- 
sult in a larger electron-capture rate during the collapse 
and a smaller trapped lepton fraction at the core bounce 
(Brucnn 1989). Throughout the birth stage of neutron 
stars including both the dynamical phase and the qua- 
sistatic phase, deleptonization will proceed faster for TMS 
than for TMI and electron-type lepton number flux will 
be larger for TMS. The work in this direction to simulate 
also the dynamical phase is in progress. 

Further studies are required to predict precisely the 
profile of the neutrino burst and to extract the information 
on the symmetry energy from the observational signals 
of the next supernova explosion in our Galaxy or close- 
by. The influence of the thermodynamical properties of 
EOS such as the symmetry energy should be carefully ex- 
amined besides the one by the stiffness, which has been 
mainly focused. The symmetry energy is one of the most 
interesting keys to clarify the current problems in nuclear 
physics and astrophysics. It would be interesting to apply 
the many body framework such as the RMF theory, in 
which the symmetry energy is being checked by the prop- 
erties of unstable nuclei, to solve consistently the current 
issues on the cooling of the neutron stars and the super- 
nova explosion together with the ones in nuclear physics. 
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Figure Captions 

Fig. 1. The energy per baryon of nuclear matter in the RMF 
theory is shown as a function of the baryon density for the cases 
with TM1 and TMS. The solid curves represent the results of 
symmetric nuclear matter and pure neutron matter with TM1. 
The dashed curve represents the result of pure neutron matter 
with TMS. The result of symmetric nuclear matter with TMS 
is the same as the one with TM1. 

Fig. 2. a The proton fraction of neutron star matter is shown 
by the solid curve for the case with TM1 as a function of the 
baryon mass density. The dashed curve represents the case with 
TMS. b The gravitational masses of neutron stars with the 
EOS of neutron star matter in the RMF theory with TM1 
(solid) and TMS (dashed) as functions of the central baryon 
mass density. 

Fig. 3. The radial profiles of the matter entropy (S) at 
t — Osec (initial models) and at t = 15 sec. Solid lines are the 
model with TM1 and dashed lines are the model with TMS. 
The abscissa is the enclosed baryon mass (Mb) in units of the 
solar mass (Mq). 

Fig. 4. The radial profiles of the temperature (T) at t = sec 
(initial models) and at t = 15 sec in the same notation as Fig. ^[ 

Fig. 5. The radial profiles of the density (pa) at t = Osec 
(initial models) and at t = 15 sec in the same notation as Fig. |§| 
The differences in density profiles between the two models are 
very small. 

Fig. 6. The radial profiles of the electron fraction (Y c ) at 
t — sec (initial models) and at t = 15 sec in the same no- 
tation as Fig. |ij 

Fig. 7. The radial profiles of the lepton fraction (Yl) at 
t — sec (initial models) and at t = 15 sec in the same no- 
tation as Fig. |§| 

Fig. 8. The radial profiles of the chemical potential of v c 
in /3-equilibrium (p„) at t = Osec (initial models) and at 
t — 15 sec in the same notation as Fig. ^. 

Fig. 9. The time profiles of the net flux of the electron type 
lepton number are shown for the cases with TM1 and TMS by 
the solid curve and the dashed curve, respectively. 

Fig. 10. The time profiles of the mean energy of i> c are shown 
for the cases with TM1 and TMS in the same notation as 
Fig. | 



Fig. 11. The time profiles of the luminosity of v c are shown for 
the cases with TM1 and TMS in the same notation as Fig. ^. 
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